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We extend our previous theoretical study of Mn-based orthorhombic metamagnets to those that 
possess large nearest neighbour Mn-Mn separations (di > 3.22 A). Based on our calculations, we 
design and synthesize a series of alloys, CoMnGei-^P^, to experimentally demonstrate the validity 
of the model. Unusually, we predict and prepare several metamagnets from two ferromagnetic end- 
members, thus demonstrating a new example of how to vary crystal structure, within the Pnma 
symmetry group, to provide highly tunable metamagnetism. 



PACS numbers: 75.30.Sg,75.30.Kz, 75.80.+q, 75.30.Et 
I. INTRODUCTION 

Manganese-based orthorhombic (Pnma) binary and 
ternary alloys are of fundamental research interest as 
they often exhibit complex, non-collinear magnetic struc- 
tures that can be tuned by temperature, pressure and 
applied magnetic field. Some well-known examples arc 
the set of fan, helical and cycloidal spin structures found 
in the (H, T) phase diagram of MnpGHS the cycloidal an- 
tiferromag netism (AFM) of IrMnSP and helical metam- 
agnetism of CoMnSP^. Various theoretical explanations 
have previously been suggested to describe the mecha- 
nisms responsible for the noncollinear magnetism of such 
materials. Some refer to com peti ng symmetric and asym- 
metric exchange interactions 2 e ; others to conduction- 
mediated indirect (RKKY) exchange^. Another poten- 
tial cause put forward is band cros sing and appropriate 
Fermi surface topology (nesting/^. 

One of the most feature-rich materials of this kind is 
CoMnSi on account of its pronounced magnetic field- 
induced tricritical metamagnetism and associated neg- 
ative magnetocaloric effect (MCE^P. Our recent high- 
resolution neutron diffraction (HRPD) studjE^ uncovered 
a giant magneto-elastic coupling within the antiferromag- 
netic ground state of this system. It occurs as a change of 
up to 2% in nearest-neighbour Mn-Mn separations d\ , e?2 
on heating. The field-induced tricriticality of this system 
can thus be understood as the result of tuning the meta- 
magnetic critical temperature with an applied magnetic 
field to the point at which it coincides with this native 
giant magneto-elasticity. 

Using Density Functional Theory (DFT) we recently 
examined the importance of the same Mn-Mn separa- 
tions in determining the occurrence of different mag- 
netic groundstates across several Mn-based orthorhom- 
bic (Pnma) systems P^l. By applying hydrostatic expan- 
sion and compression to a prototype model MnP al- 
loy, we found a stability criterion for the appearance of 
an AFM groundstate, rather than the usual FM state 
seen in MnP. This direct relation between Mn-Mn sep- 
aration and magnetic groundstate can explain the en- 
ergetic proximity of FM and AFM states in materi- 
als such as CoMnSPS, MnA Sl _ x Pji2l, (Fei_ x Co x )MnF<S] 



and NiMnGei- x SiJ^ where the nearest-neighbour Mn- 
Mn distances are close to a critical separation of 2.95 A< 
di <3.05 A. 

Although the latter model was computed in a large in- 
terval of 2.5 A< d 1 <3.22 A to cover many of the relevant 
compositions in the literature, it may lead to a misinter- 
pretation of a technologically relevant alloy with a larger 
di = 3.4 A spacing: CoMnGe. CoMnGe is a collinear 
ferromagnet with a tendency to form a metastable hexag- 
onal structure upon rapid cooling^. In this Article, we 
extend our previous theoretical analysis^ towards larger 
Mn-Mn separations to explain the re-appearance of fer- 
romagnetism in alloys with d\ >3.37 A. The importance 
of the correct theoretical description of latter composi- 
tion is due its large magnetocaloric effect around room 
temperature^^. 

In this Article we first show the striking re-appearance 
of a FM groundstate at large interatomic Mn separations 
in Mn-based Pnma alloys where d\ >3.37A, thereby ac- 
counting for the magnetic properties of CoMnGe. Sec- 
ondly, and significantly, we have designed a new alloy se- 
ries, CoMn(P,Ge) based on our extended model in order 
to test and demonstrate its validity and in particular the 
dominance of the (di) Mn-Mn separation in determining 
the magnetic groundstate of the series of Mn-containing 
Pnma alloys. We show that metamagnetism can be de- 
rived, unusually, by inter-doping two ferromagnetic end- 
compositions in order to bring d\ to the critical regime 
where antiferromagnetism and ferromagnetism are simi- 
lar in energy. 

The remainder of the Article is organised as follows: 
first in Sec. [TTJ the theoretical results calculated by ap- 
plying DFT to the prototype MnP alloy are given. Based 
on this model, we present the structural and metamag- 
netic properties of purposefully designed CoMnGei-^P.,; 
alloys in Sec. |III| Finally, a summary is made and con- 
clusions are drawn in Sec. IIVI 



II. THEORETICAL 

In our previous work we considered what we term the 
"prototype" binary MnP (Pnma) alloy and calculated 
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the effect of isotropic lattice expansion and compression 
on hypothetical non-magnetic (NM), ferromagnetic (FM) 
and antiferromagnetic (AFM) states by using the gen- 
eral gradient aproximation method (GGA-DFT) imple- 
mented in VASP Kresse and Furthmuller 19 . We found 
the critical lattice parameters where a crossover from one 
magnetic state to another can occur^S In that study, a 
single unit cell consisting of 8 atoms (4 Mn and 4 P) was 
used, which allowed three different collinear antiferro- 
magnetic configurations (AFM1, AFM2 and AFM3) and 
a collinear ferromagnetic (FM) one to be constructed. In 
the interval of 2.5 A< d x <3.22 A we were able to predict 
a transition in the zero temperature magnetic structure 
from NM to FM, and finally to AFM as a function of 
expanding lattice parameters. A detailed description of 
the DFT calculations is given in that wortP. Here we 
extend our simple binary model to g?i>3.22 A values by 
further hydrostatic expansion in order to interpret ferro- 
magnetism in CoMnGe where d\ = 3.4 A in the current 
model. 

Fig. ^ shows the difference in energy between the pos- 
sible collinear AFM and FM magnetic states (AExot — 
Eafm — Efm) as a function of Mn-Mn separation. Using 
this comparison scheme a non-FM state becomes most 
favourable when it has the most negative value of AExot- 
On the left hand side of Fig. [T] the large compression 
causes a strong overlap of d-orbitals, and the broad d-d 
hybrid bands thus formed cannot support spontaneous 
magnetisation. In agreement with experimental findings, 
the FM state is the groundstate of MnP and is stable 
for intermediate deformations of the lattice. On fur- 
ther expansion of the lattice parameters, first the AFM1 
type magnetic structure (at around c£i~2.97 A) and then 
the AFM3-type ordering (at d±~3.1 A) become energet- 
ically favorable. However the extended study presented 
here shows that AFM3-type ordering ceases to be the 
most stable magnetic state for large lattice expansion, 
and eventually the collinear FM state is once again the 
groundstate for di>3.37 A. 

In order to experimentally prove the validity of the 
model, we carefully selected two collinear FM compo- 
sitions with lattice parameters from the different FM 
regions of the stability plot in Figure [T] Our choices 
were CoMnP from the FM(1) regime (di~2.95A>22l) and 
CoMnGe (^^3.4^) from the larger Mn-Mn separation 
(FM2) zone. Our hypothesis is that progressive substitu- 
tion of one p element for another in CoMnGei_a;P x , with- 
out changing the 3d-element concentration, can cause the 
appearance of metamagnetism in particular compositions 
of the series. This replacement of large germanium atoms 
by the much smaller phosphor atoms should result in a 
shrinkage of lattice. From Figure [I] we expect that the 
decreasing di separation will lead to the destabilization 
of FM state in competition with the AFM one at a cer- 
tain P/Ge ratio. It should then be possible to manipulate 
the magnetic state of the energetically metastable alloys 
by changing temperature or applied magnetic field. 

In order to identify the key factors that can lead to 
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Figure 1: (Color online) Stability of possible collinear mag- 
netic structures, relative to ferromagnetism, within a single 
unit cell of MnP as a function of di Mn-Mn separation. AFM 
configurations become stable where AEafm-fm < 0. The 
vertical dashed line represents the experimental (strain-free, 
e = 0) lattice of MnP. We see that the FM state is first desta- 
bilised by lattice expansion, and then becomes stable again 
at large d\ values: FM(2). 



the substantial changes in magnetic groundstate, we first 
calculate the electronic band structure and magnetic mo- 
ment of CoMnGei^Pz alloys with x=0, 0.25, 0.5, 0.75 
and 1 in a collinear FM state. The calculated magnetic 
moments are given in Table [I] The magnetic moments on 
Mn (2.68^b) and Co (0.28/x B ) sites in the CoMnP alloy 
agrees well with previous calculations based on the KKR 
method with coherent potential approximation (CPA) 
by Zach and co-workers2B. Furthermore, the partial re- 
placement of P by Ge results in a progressive increase 
of magnetic moment on both 3d elements, leading to an 
increased Mxot of up to 3.58 /ib for CoMnGe. A small 
negative moment induced on the p-block elements is also 
observed. 

The FM total density of electronic states (DoS) of the 
two end compositions, CoMnP and CoMnGe, are plot- 
ted, together with CoMnP .5Ge .5 in Fig. [2] Although 
the value of total density of states at the Fermi level, 
NTot(Ev) — N^(Ef) + N-ff(Ef) exhibits a large change 
with composition, each total DoS possesses the same fea- 
tures over a large extent of energy range. The main dif- 
ference is the location of these features, and in particular 
the location of a pseudogap-like feature in the DoS. Us- 
ing CoMnP as reference, if the energy scale of the minor- 
ity DoS is shifted by about +0.25 eV for CoMnGe and 
about -0.3 eV for CoMnP .5Ge .5, not only would the 
pseudogap fall at Ef but most of the exchange-split DoS 
peaks of Mn and Co would line up at around the same 
position. The large NTot{Ep) in CoMnP .5Ge .5 (Table 
[l| suggests the instability of the collinear FM state in 
this composition. A possible scenario that stabilizes the 



3 




CoMnGe P 

0.75 0.25 



E-E F (eV) 



Figure 2: (Color online) Collinear FM total density of states 
for CoMnGei-^Pa; with x=0, 0.5 and 1. The Fermi energy 
falls into the hybridization gap for both CoMnP and CoMnGe, 
but shifted for CoMnPo.sGeo.s resulting in a large NTot(Ef). 
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Table I: Calculated magnetic moments (/xs) and NTot(Ef) 
(states/eV/f.u.) for CoMnGei-^P^. 



noncollinear state through the formation of hybridization 
gap at the Fermi energy is described by Lizarraga et al.^. 
We recently showed the relevance of this mechanism in a 
noncollinear DFT study on the metamagnet CoMnSf^l 
In the following sections, we are going to demonstrate 
the validity of our theoretical prediction of metamag- 
netism in these Mn-based Pnma alloys through magnetic 
and structural results on an experimentally synthesized 
CoMnGei-^P^ series of alloys. 



III. EXPERIMENTAL 

A. Experimental details 

Samples of CoMnGei-^Pj; with x = 
0.25,0.4,0.5,0.55,0.6 and 0.75 were prepared in a 
quartz nozzle by an induction melting technique, using 
Co 2 P (99.9%) and Mn 3 P 2 (99.9%) master alloys mixed 
together with high purity Co (99.97%), Mn (99.99%) 
and Ge (99.9999%) elements in the required propor- 
tions. The alloys were cast into a copper mold under 
protective Ar atmosphere. The ingots were then sealed 




Figure 3: Representative SEM micrograph and corresponding 
elemental mapping from EDX of CoMnGeo.75Po.25- 



in quartz tube under protective He atmosphere and a 
homogenization at 1000°C for 24 hours and annealing 
treatment at 800°C for 72 hours then followed. The 
samples thus obtained were crushed into fine powder in 
order to determine their crystal structure using X-ray 
diffraction (XRD) with Cu Ka radiation. Structural 
(Rietveld) refinement of the data was carried out the 
using the Fullprof ?? program. A microstructural and 
compositional analysis was carried out using a Carl Zeiss 
1540EsB scanning electron microscope (SEM). Magnetic 
properties of the samples were studied in a Quantum 
Design MPMS system. 



B. Crystal Structure 

Both CoMnGe and CoMnP alloys crystallize in the or- 
thorhombic (Pnma, 62) structure in which the elements 
occupy general 4c (x,^,z) crystallographic positions. 
XRD analysis of the CoMn(Ge,P) samples revealed the 
formation of this orthorhombic structure in all composi- 
tions. Furthermore, in the samples with x = 0.25, 0.4, 0.5 
and 0.75, extra reflections in the diffraction pattern also 
appear that can be ascribed to the hexagonal Ni2ln-type 
(P63/mmc, 194) lattice structure. The appearance of 
the higher symmetry hexagonal phase is often observed 
in similar alloy systems because the orthorhombic struc- 
ture can be regarded as a distortion from this hexagonal 
structure and the two structures can be interrelated as 
follows: b or tho = a he x and c ortho — %/3 x a hex . 

The importance of this latter correlation has 
been exploited in several Mn-based Pnma sys- 
tems. In the pseudo-binary Mn!_ x Fe x As alloys the 
sharp, first order type magnetostructural transforma- 
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Figure 4: (Color online) Volume fraction of the orthorhombic 
CoMnGei-xPj and hexagonal CoMnGe phases as a function 
of nominal composition (top) and lattice parameters of the 
main orthorhombic phase as a function of corrected phosphor 
concentration (bottom). The lattice parameters of CoMnSi, 
a metamagnet, are added (star symbols) for comparison. 



ference between the two structures. These findings are 
direct evidence for a compositional phase separation of 
the quaternary alloy rather than for a second-order type 
transformation of the single quaternary composition with 
temperature. 

The results of the quantitative Rietveld analysis in Fig. 
[4] (top) reveal the formation of single phase orthorhombic 
structures in the alloys rich in P (ir>0.5). Fig. [4] (bot- 
tom) summarizes the lattice parameter of the orthorhom- 
bic structure as a function P content. (The composition 
values (x 1 ) given in the lower figure are corrected based 
on the quantitative analysis.) The b and c lattice con- 
stants show decrease continuously with Ge addition while 
the a parameter stays nearly constant and until the lat- 
tice parameters with x~0.5 become close to those of the 
metamagnet CoMnSi. 

In the next section, we will demonstrate the occurrence 
of metamagnetism in these Mn-based alloys that have 
appropriately designed lattice parameters. 



C. Magnetic properties 



tion (orthorhombic<^=> hexagonal p21 is also accompanied 
by a "colossal" MCEP. The magnitude of the use- 
ful magnetic entropy change is, however, now strongly 
contested^. A similar magnetostructural transition in 
CoMnGe-based ternary compositions was reported by 
Kanomata and co-workers who observed a large, ~5.3%, 
volume change^2ZI. Theoretical calculations also re- 
vealed that the Co vacancy-induced phase transforma- 
tion is due to a high moment to low moment magnetic 
transition accompanied by a large magnetovolume effect 
originating from the change of the coupling distance be- 
tween the principal magnetic atoma ^. A s an extension of 
this study, both HamerSSbnd TrunjpESl recently demon- 
strated that the transitions can be fine-tuned by pseudo- 
ternary additives (Sn, B or Cr elements) around room 
temperature for an enhanced magnetocaloric effect. The 
appearance of the hexagonal structure at room temper- 
ature in our samples can therefore be understood as a 
first order transformation from the low temperature or- 
thorhombic structure to a high temperature hexagonal 
one. 

In order to clarify the composition of the hexagonal 
phase, a representative secondary electron (SE) SEM mi- 
crograph was taken from the CoMnGeo.75Po.25 alloy and 
is shown in Fig. [3^,. The two different structural regions 
are identified in accordance with the XRD results. EDX 
elemental mapping revealed that the two phases have 
substantially different atomic compositions. Our inves- 
tigations show that the orthorhombic phase is poor in 
Ge (and consequently enriched in P) whereas the second 
phase is enriched in Ge and is therefore close in composi- 
tion to stoichiometric CoMnGe. An atomic composition 
profile taken along the direction of the arrow indicated 
in Fig. [5J) shows quantitatively the compositional dif- 



As we demonstrated in Sec. |III B| the lattice param- 
eters of the CoMnPi_ 2; Ge. r alloy can be tuned towards 
those of the CoMnSi metamagnet. In the present sec- 
tion, we will show that this structural engineering allows 
us to prepare metamagnetic quaternaries, even though 
the end alloys (CoMnP and CoMnGe) are ordinary fer- 
romagnets. Based on the a, b and c lattice parameters of 
the CoMnGei-^Pa; alloys, one would expect progressive 
change from collinear ferromagnetism to non-collinear 
anti-ferromagnetism with the replacement of P by Ge in 
the system. At around x~0.5 the structure should resem- 
ble that of the CoMnSi metamagnet and once the compo- 
sition is on the P-rich end of the series, ferromagnetism 
should re-appear. Indeed, we have found this scenario to 
be fulfilled; however the structural phase separation into 
a ternary CoMnGe and a residual CoMn(Ge,P) phase in 
compositions that are rich in Ge hinders a simple inter- 
pretation. 

The composition dependence of magnetisation loops 
taken at 10 K is shown in Fig. [5] Metamagnetism is most 
visible in the single phase compositions with x=0.55 and 
0.6. In these samples, the initial magnetisation varies al- 
most linearly with applied field below the inflection point 
that occurs at a critical field value. Although this critical 
behaviour is also visible in all compositions with x < 0.6, 
the increasing volume of ferromagnetic CoMnGe phase 
at low x suppresses the sharpness of the upturn in mag- 
netisation. A distinct shift of the peak of seen in 
Fig. [5]d shows that the critical field decreases sharply as 
more Ge is replaced by P (towards small x). 

The composition dependence of the critical field in a 
field of 1 Tesla is shown in Fig. [6] and is consistent with 
the magnetisation loops collected at 10 K. The lower the 
critical field that is observed at 10 K (Fig. |5|, the lower 
the critical temperature is at which the 1 T applied mag- 
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Figure 5: (Color online) Magnetisation loops of 

CoMnGei-^ at 10K. 



netic field is sufficiently large to bring the sample to the 
high magnetic state. A strong dependence of the meta- 
magnetic transition temperature on crystal structure is 
evident from this Figure, when compared with Fig. [4] 
Although samples with x < 0.5 contain a minor sec- 
ond phase, they all show a characteristic increase in the 
magnetisation as a function of temperature, indicative of 
native metamagnetism. In the high Ge (low x) side of 
the series, with 2=0.25, although the substantial ferro- 
magnetic volume (from the CoMnGe phase) largely sup- 
presses the metamagnetic transition, it still appears at 
around room temperature (RT). 

As the samples become richer in phosphor, the criti- 
cal temperature decreases sharply from ^300K to ~80K 
and it eventually disappears for 2=0.75. The highest 
magnetisation values in a 1 T applied field show an in- 
crease with x (except for 2=0.25) as the result of a bal- 
ance of several effects. Firstly, it is easily foreseeable that 
at lower transition temperatures ferromagnetic configu- 
rations will exhibit larger overall net moments compared 
to those of the higher temperature ones. Secondly, the 
presence of the hexagonal CoMnGe lifts the low temper- 
ature "baseline" of magnetisation in Fig. [6j 



IV. SUMMARY AND CONCLUSIONS 

Using DFT calculations based on a "prototype" bi- 
nary MnP composition, we have investigated the occur- 
rence of AFM and FM states in Mn-based orthorhombic 
(Pnma, 62) alloys. As the result of isotropic expansion, 
the FM(1) ground state is no longer stable but instead 
AFM coupling of the spins on the Mn atoms is predicted 
above d\ >2.95 In this work, we have extended our 
theoretical investigation to higher hydrostatic expansions 
and found the re-occurrence of ferromagnetism (FM2) 
at large Mn-Mn separations over di >3.37 A that also 
explains the collinear ferromagnetism in CoMnGe with 
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Figure 6: (Color online) Iso-field magnetisation of 
CoMnGei-^P^ as a function of temperature, measured in a 
1 Tesla applied field. We observe metamagnetism in almost 
all samples, and a sharp change in the metamagnetic critical 
temperature with composition in the range £=0.4 to 0.6. 



di = 3.4A. 

Based solely on the deduced magnetic stability plot 
(Fig. [IJ, we designed a series of pseudo-ternary CoMn- 
based alloys in order to experimentally prove the valid- 
ity of our theoretical concept. Taking two collinear FM 
ternaries: one, CoMnP with a low d\ from the FM1 re- 
gion ; the other, CoMnGe with a high d\ from the FM2 
region we attempted to drive the alloy magnetism to- 
wards the metamagnetic/ AFM zone by careful structural 
design. 

The experimental investigation of CoMnGei_ ;r P. E has 
indeed revealed an AFM ground state for compositions 
x w 0.5. The appearance of a magnetic field- and 
temperature-dependent metamagnetic transition in sev- 
eral samples also suggests the existence of complex non- 
collinear spin structure in most of them , and in par- 
ticular in the range 2=0.4 to 0.6. The large predicted 
NTot{Et) for x = 0.5 (in Sec. [TTJ) in a hypothetical 
FM state is because of a shift in the energy of the 
hybridisation-derived pseudogap as the lattice parame- 
ters expand upon Ge substitution for phosphor. 

Although the complex magnetic spin structure of these 
new samples is to be determined, the system is an ex- 
ample of the stabilization of non-collinear magnetism 
through the formation of hybridization gap at the Fermi 
energy as described by Lizarraga et al.^ and as recently 
found in the ternary CoMnSp^l. The Mn-containing 
Pnma structure is extremely versatile with regard to 
elemental substitution. The above demonstration of a 
structurally-directed tuning of magnetic properties there- 
fore provides a potential direction for future tailoring of 
metamagnetic phase transitions towards their use in ap- 
plications such as those that rely on the magnetocaloric 
effect. 
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